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\0 : ABSTRACT 

We have investigated the spatial extent and structure of the Draco dwarf spheroidal galaxy 
^ ■ using deep wide-field multicolor CCD photometry from the Sloan Digital Sky Survey (SDSS). 

CD ! Our study covers an area of 27 square degrees around the center of the Draco dwarf and reaches 

■ two magnitudes below the level of Draco's horizontal branch. The SDSS photometry allows very 
effective filtering in color-magnitude space. With such filtering the density of the foreground of 
Galactic field stars is decreased by more than an order of magnitude and the stellar population of 
the Draco dwarf galaxy thus stands out with much higher contrast than in former investigations. 

■ We show that the spatial distribution of Draco's red giants, red horizontal branch stars and 
subgiants down to i*=21.7 mag does not provide evidence for the existence of tidally induced 
tails or a halo of unbound stars. The projected surface density of the dwarf galaxy is flattened 
with a nearly constant cUipticity of 0.29 ± 0.02 at position angle 88 ± 3°. The radial profile can 
be fit by King models as well as by a generalized exponential. Using the empirical King (1962) 

^ I profile the core radius and the limiting (or tidal) radius along the major axis are Tc = 7'.7 and 

i rt — 40.'1, respectively; the latter means that the size of the Draco dwarf galaxy is 40% larger 

than previously measured. Fitting the profile of King's (1966) theoretical models yields a still 
larger limiting radius of rt = 49f5. There is no clear indication of a tail-like extension of the Draco 
population beyond this radius. A break in the radial surface density profile which might indicate 
a halo of extratidal stars is also not found in our Draco data. We conclude that down to the 
above magnitude limit tidal effects can only exist at a level of 10"'^ of the central surface density 
of Draco or below. The regular structure of Draco found from the new data argues against it 
being a portion of an unbound tidal stream and lends support to the assumption of dynamical 
equilibrium which is the basis for mass estimation. The changes in the values for the core radius 
and limiting radius imply that the total mass of Draco is higher by more than a factor of two. 
Using a King (1966) spherical model of equivalent size as a reference and adopting a line-of-sight 
velocity dispersion of cither 10.7 kms~^or 8.5 kms~^(Armandroff et al. 1995), we derive estimates 
of the total mass within radius of 3.5 ± 0.7 x 10^ Mq and 2.2 ± 0.5 x 10^ Mq, respectively. 
From the combined i band flux of all possible Draco members that lie within major axis radius 
Tf we determine the total luminosity of the Draco dwarf galaxy as {L/LQ)i = 2.4 ± 0.5 x 10^. 
This includes corrections for the flux of the foreground stars and the unseen fainter part of the 
Draco population. We thus obtain overall mass-to-light ratios M/Li of 146 ± 42 or 92 ± 28 in 
solar units. In summary, our results strengthen the case for a strongly dark matter dominated, 
bound stellar system. 

Subject headings: galaxies: individual (Draco) — galaxies: dwarf — galaxies: structure — galaxies: 
Local Group 
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1. Introduction 

The Milky Way is surrounded by a number 
of small faint companion galaxies that are gravi- 
tationally bound to it. Nine such Milky Way 
dwarf satellites at galactocentric distances be- 
tween 16 kpc and 250 kpc are presently known. 
The dwarf spheroidal galaxy in Draco (hereafter 
called Draco dSph or simply Draco), located at 
a distance of about 71 kpc (see § 5.1), is one of 
the two least luminous of these companions. Esti- 
mates of the dynamical masses of the Milky Way 
dwarf satellites based on their spatial extent and 
internal velocity dispersion lead to the conclu- 
sion that most of tlicnn have high mass-to-light 
ratios (30 - 100) and hence must contain large 
amounts of dark matter (see Mateo 1998 and ref- 
erences therein). The parameters of the Draco 
dSph imply a particularly high mass-to-light ra- 
tio of M/L ~ 60 - 90 (Armandroff et al. 1995). 
Thus Draco must be completely dominated by 
dark matter if the underlying assumption of dy- 
namical equilibrium holds. 

However, even with a large fraction of dark 
matter the total masses of dwarf spheroidal galax- 
ies like Draco are at most a few 10^ Mq while their 
spatial extent is several hundred parsecs. There- 
fore these systems are susceptible to perturbations 
by the tidal forces of the Galactic potential if their 
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orbits are sufficiently eccentric to bring them to 
galactocentric distances of about 20 kpc or closer. 
As the result of such perturbations, a dwarf galaxy 
may lose significant amounts of mass and eventu- 
ally be disrupted. This has a number of important 
implications: 

1. Tidal mass loss and disruption of dwarf galax- 
ies can shed light on the formation of the stellar 
Galactic halo and its spatial and kinematic struc- 
ture (Johnston et al. 1996, Johnston 1998, Helmi 
k White 1999). 

2. The spatial distribution of tidal debris can 

provide constraints on the space motions and or- 
bits of the companions (e.g., Odenkirchen et al. 
2001). These would, e.g., be useful for testing 
the hypothesis that the Milky Way satellites are 
grouped in certain streams (e.g., Lynden-Bell & 
Lynden-Bell 1995). Traces of the orbits of the 
dwarf galaxies would also provide constraints on 
the gravitational potential of the Galactic halo 
(e.g., Kuhn 1993, Johnston et al. 1999a). 

3. If a dwarf galaxy is tidally perturbed it may 
be incorrect to estimate its mass via the assump- 
tion of dynamical equilibrium. Kuhn & Miller 
(1989) have argued that tidal heating can inflate 
the velocity dispersion of a dwarf galaxy and hence 
lead to overestimated masses and mass-to-light ra- 
tios. Numerical simulations by Klcssen & Kroupa 
(1998) have shown that in the extreme case of 
complete disruption the observed properties of lo- 
cal dwarf galaxies may under certain conditions 
be explicable without any dark matter. This point 
however is controversial (see Olszewski 1998 and 
references therein). Piatck & Pryor (1995) found 
that even when tides affect the structure and the 
kinematics of a dwarf spheroidal, they fail to in- 
flate the central mass-to-light ratio. 

An impressive example of tidal interaction is 
given by the Sagittarius dwarf spheroidal, the clos- 
est companion of the Milky Way which is observed 
while being disrupted in a close encounter with 
the Galaxy (Ibata et al. 1997). Sagittarius has 
a very elongated and distorted shape and part 
of its stellar content is stretched out in a tidal 
stream along its orbit (Mateo et al. 1998, Majew- 
ski et al. 1999, Johnston et al. 1999b, Ibata et al. 
2000). This shows that tidal effects can indeed 
play an important role. It is therefore desirable to 
find out if the other (more distant) dwarf satellites 



2 



also reveal signs of tidal mass loss and disruption. 
Majewski et al. (2000) showed that the giants of 
the Carina dwarf spheroidal have an extended dis- 
tribution that rc'ac;lK's beyond the estimated tidal 
limit and produces a break in the projected radial 
density profile. Martinez-Delgado et al. (2001) re- 
cently reported the detection of a tidal extension 
for the dwarf spheroidal in Ursa Minor. 

Prom recent studies of Draco there seemed to 
be evidence that this dwarf also has an extrati- 
dal stellar component. Irwin & Hatzidimitriou 
(1995; hereafter IH95) investigated the spatial 
structure of all known dwarfs except Sagittarius 
by means of star counts from digitized photo- 
graphic plates. They found for Draco, as well 
as for most of the other dwarfs, significant over- 
densities of stars outside the limiting radius of 
their best-fit King profile. Smith, Kuhn & Haw- 
ley (1997, hereafter SKH97) reported preliminary 
results of CCD photometry in a number of fields 
along the major and minor axis of the Draco dSph 
and claimed to see extratidal members out to 3° 
angular distance from the center of Draco in east- 
ern direction. Piatek et al. (2001) investigated 
Draco by analysing color-magnitude diagrams of 
nine CCD fields at different positions and angular 
distances from its center. They found evidence for 
Draco stars in those fields which overlap with the 
standard limiting radius of 28.'3 given by IH95 and 
reach somewhat beyond this limit. Furthermore, 
they argued that there might be weak evidence 
for Draco stars also in their more distant fields at 
about 1.3 degrees from the center of Draco. 

The Sloan Digital Sky Survey (SDSS; see York 
et al. 2000, Gunn et al. 1998, Fukugita et al. 1996) 
has recently enabled us to investigate the Draco 
dSph through multi-color CCD photometry in a 
large contiguous field. Hereby, the distribution of 
its stellar content has been mapped more compre- 
hensively and in more detail. As we show in this 
paper, the SDSS data do not confirm the conclu- 
sions from previous studies. 

The paper is organized as follows: § 2 provides 
details of the SDSS data set and describes our 
method for the photometric selection of Draco 
members. In § 3 we present surface density maps 
for two samples of member candidates and inves- 
tigate them for signs of tail- like extensions. In § 4 



wc discuss fits of ellipsoidal models to the observed 
distribution of stars and determine the basic spa- 
tial parameters of the galaxy. § 5 presents revised 
estimates of the total mass, luminosity and mass- 
to-light ratio. In § 6 we summarize the results and 
conclusions. 



2. Photometric selection of Draco candi- 
dates 

SDSS is a multipurpose project which is de- 
signed to provide homogeneous, deep CCD imag- 
ing in five passbands {u,g,r,i,z) with contiguous 
coverage of 10^ square degrees in the Northern 
Galactic Cap. The SDSS data are thus ideally 
suited for studies of Galactic structure and local 
stellar systems (see, e.g., Yanny et al. 2000, Ivezic 
et al. 2000, Odcnkirchen et al. 2001, Chen et al. 
2001) even though the primary goal of SDSS lies 
in extragalactic and cosmological research. 

2.1. The data set 

The field of Draco was observed in four neigh- 
boring great circle scans made with the SDSS 
CCD mosaic camera on April 3/4, and April 5/6, 
2000 (SDSS runs 1336, 1339, 1356, and 1359). The 
observations cover a four degree wide band across 
the Draco dSph galaxy and its surroundings and 
reach down to a limiting magnitude of r* w 23.0. 
Our study focusses on the declination range of 
55° < ^ < 61° and on objects recognized as point 
sources (based on their images being consistent 
with the local point spread function (PSF)). The 
so-defined sample of about 290,000 objects has 
been investigated using PSF photometry from the 
standard SDSS data reduction pipeline (Lupton 
et al. 2001). The typical photometric accuracies 
in the data set (median internal errors) are better 
than 0.03 mag in g, r, and i for sources brighter 
than 20.0 mag in these filter, and reach the level 
of 0.1 mag for source brightnesses of u* = 21.0, 
g* = 22.3, r* = 22.0, i* = 21.5, and z* = 20.0. 
The data are preliminary in the sense that the 
calibration of the instrumental system onto to 
the nominal SDSS photometric system is still in 
progress. Therefore, we use asterisks to designate 
the currently used system of magnitudes (York et 
al. 2000). The differences that may occur with 
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Table 1 

Basic parameters of the Draco dsph galaxy 



Parameter 


Value 


Reference 


a, 5 (J2000) 


17''20™13?2, +57°54'54" 




l,b 


86.4°, +34.7° 




d 


71 kpc 


Stetson (1979) + Carney et al. (1992) 


Vr 


-293.3 kms-l 


ArmandrofT ct al. (1995) 


<^los 


10.7 kms-i 


Armandroff et al. (1995) 


[Fe/H] 


-2.0 dex 


Lehnert et al. (1992) 



respect to a later final calibration will however be 
small (not larger than a few percent) and will thus 
not affect the results of this study in a significant 
way. 

The IR dust emission maps of Schlcgcl ct al. 
(1998) show that some of the regions around the 
Draco dSph galaxy have non-negligible and spa- 
tially varying extinction. This is potentially dan- 
gerous for our goal of detecting extended struc- 
tures in the Draco population because it creates 
spatially varying shifts in the colors of the stars. 
The SDSS data include for each object and each 
filter an individual extinction estimate which is de- 
rived from the reddening proposed by the Schlcgcl 
et al. (1998) maps. Depending on the filter and 
the position of the object, the estimated extinc- 
tion within our sample varies between 0.03 and 
0.3 mag in g, r and i, with local peaks of 0.4 
in u. To remove reddening effects, the estimated 
extinction values have been subtracted from the 
measured magnitudes prior to any further anal- 
ysis. We show the estimated extinction, Ag, in 
Fig. 1 as a function of position on the sky. For the 
main body of Draco, Ag is, on average 0.10 mag, 
while west of the galaxy and towards the southeast 
and the east there arc regions where Ag increases 
above 0.16 mag or 0.2 mag. The latter corresponds 
to a shift of about 0.05 mag in the color index g—i. 

The crucial point in constructing a high- 
contrast map of a stellar system like Draco is 
efficient discrimination between the tracers of this 
system and the foreground field stars. Therefore, 
specific photometric selection criteria had to be 
worked out, using the central part of Draco as 
a training set. In order to determine the pho- 



tometric characteristics of the Draco population 
in the SDSS system we analysed the stars in an 
ellipse of semi-major axis 9' around the nominal 
center. These stars shown in Fig. 2. The mem- 
bers of the Draco dSph have colors in the range 
of -0.5 < g* -r* < 1.0 and -0.5 < r* - i* < 0.5 
(see panels (a) and (b) of Fig. 2). The uncertain- 
ties of the g* , r* and i* magnitudes of these stars 
are approximately equal while the uncertainties of 
their u* and z* magnitudes are at least a factor 
of two higher. Thus the g*, r*, and i* magni- 
tudes were identified as the most siiitable quanti- 
ties for photometric filtering. The Draco members 
show a strong correlation between the color indices 
g* — r* and r* — i* which can be expressed in a 
simple linear relation between these indices. This 
fact allowed us to facilitate the selection of Draco 
candidates by an orthogonal transformation from 
g* — r* , r* — i* to new indices ci, C2 as given in 
Eq.(l), similar to the indices used by Odenkirchen 
et al. (2001). 

ci = 0.927(.g*-r*) + 0.375(r*-7;*) (la) 
C2 = -0.375(5* - r*) + 0.927(r* - i*) (lb) 

Panels (c) and (d) of Fig. 2 present color- 
magnitude diagrams for the central part of Draco 
using these indices. The plane of (ci,i*) shows 
the familiar pattern of an old stellar population, 
i.e. red giant branch, subgiant branch and hori- 
zontal branch. In the plane of {c2,i*) the Draco 
stars simply scatter around C2 — 0. This scatter 
is mainly due to photometric errors. An adequate 
selection in C2 was thus obtained by deriving the 
median error of C2 from the individual er- 
ror estimates in g*, r* and i* as a function of 
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source brightness i* and by imposing the con- 
straint |c2| < 2 a 02(1*) ■ The distributions of the 
color indices u* — g* and i* — z* were found to be 
rather diffuse and turned out to be of httle use for 
the selection of Draco candidates. Thus, for the 
reason of clarity, these indices were not included 
into the selection process. 



2.2. The color-magnitude mask 

The most important means to discriminate 
Draco members from field stars is the combina- 
tion of source brightness {i* in our case) and color 
index ci. Our approach was to design an em- 
pirical filter mask in the color-magnitude plane 
of (ci,i*) with optimal adaptation to the color- 
magnitude distributions of the Draco population 
and the field stars. Fig. 3 shows the density dis- 
tribution fu of Draco stars in the plane of (ci,i*) 
(Fig. 3a) and the corresponding distribution fp 
for field stars (Fig. 3b). These distributions were 
obtained through counts on a 0.01 mag x 0.05 mag 
grid and subsequent weighted averaging within a 
radius of five grid points. To get reliable count 
rates and a truly representative sample the dis- 
tribution was derived from a relatively large 
area, i.e. an ellipse with semi-major axis 30' (el- 
lipticity 0.3, position angle 90°) around the cen- 
ter of Draco. The non-neglible, but well-defined, 
contribution from field stars in this area was sub- 
tracted by means of fp- The distribution fp was 
obtained from stars with > 2° angular distance 
from the center of Draco. It exhibits two main 
concentrations: The brighter one (max. density 
at i* < 18, ci « 0.4) consists of slightly evolved 
thick disk stars and has little photometric overlap 
with the Draco population, the fainter one (max. 
density at i* > 19, ci « 0.3) consists of halo stars 
at and somewhat below the turn-off point of the 
subdwarf main-sequence, (see Chen et al. 2001).^^ 
We assume that the Draco population and the field 
population follow the distributions Jd and fp in 
all parts of the field. This is presumably true for 
the field stars, but not necessarily for the Draco 
stars if there are population gradients. Neverthe- 
less, one must adopt homogeneity as a working 



hypothesis. In Fig. 3c we show the population 
contrast s = fo/fr as a function of position in 
the plane of (ci.i*). The maxima of s delineate 
those regions in the color-magnitude plane with 
the highest fraction of Draco stars. To construct a 
map of Draco, both high contrast s and large num- 
bers of stars are necessary. This suggests that a 
filter mask for the selection of Draco stars should 
include all points with s > sq, with sq being a 
threshold that needs to be optimized. Variations 
in So define a one-parameter family of filter masks, 
which produce samples of different size and com- 
position. Supplementary choices had to be made 
to fix the lower limit of ci and the upper limit of 
i*, which are not well constrained by the distri- 
bution s. For the color index we chose a cut-off 
at ci = 0.0 mag since there are not many Draco 
stars beyond this limit according to Fig. 2c. In i* 
wc imposed a limit of i* < 21.7 mag, because the 
data from different observing runs have somewhat 
different limiting magnitudes and completeness 
limits that begin to show up by spatial inhomo- 
geneities beyond i* = 21.7 mag. 

To estimate the optimal threshold, sq, we com- 
puted the signal-to-noise ratio (SNR) for the num- 
ber of Draco stars from basic Poisson statistics of 
number counts with background subtraction: 



SNR(so) 



N{sa) -wNp{so) 
y/N{so)+wmp{so) 



(2) 



^^A third major field component would be M-dwarfs of the 
thin and thick disk. These do not appear in the diagram 
because they have ci > 1.0 and hence are of no concern. 



Here, N is the total number of stars in the 
sample defined by sq and within a given region of 
Draco, Np the number of stars in the same sam- 
ple, but in the region where the foreground field 
population is probed, and w scales the angular ex- 
tent of these two regions. Since the main goal was 
to determine the structure of the outer part of the 
Draco dSph, the SNR was computed and maxi- 
mized for the annulus between 25' and 75' angu- 
lar distance from the center of Draco. The max- 
imum SNR was determined by scanning through 
a dense series of threshold values sq. In Fig. 3c 
we show the color-magnitude masks in three typi- 
cal cases: the mid-level contour corresponds to the 
mask that yields the highest SNR. The masks that 
correspond to the low-level and high-level contours 
of Fig. 3c yield larger and smaller samples, respec- 
tively, but both with lower SNR. 
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3. The spatial distribution of "Draco-like" 

stars 

We now discuss the distribution of Draco can- 
didate stars in the plane of the sky by means of 
two special samples. The first sample (SI) has 
been obtained with the suboptimal filter given by 
the lowest contour of Fig. 3c; the second one (S2) 
has been selected with the optimal filter given by 
the mid-level contour of Fig. 3c. SI has twice 
as many field stars and 1.4 times as many Draco 
stars as S2. The filters reduce the surface density 
of field stars from 2.57 arcmin"^ (unfiltered case) 
to 0.15 arcmin-2 (SI) and 0.074 arcmin'^ (S2), 
i.e. by factors of 17 and 35, respectively. In turn, 
the contrast between the Draco population and 
the residual field star population is strongly en- 
hanced. Without any filtering the central stellar 
density contrast would be 5.7 (see Fig. 1), whereas 
for SI and S2 it is 68.7 and 96.9, respectively. 
Our samples thus reveal the spatial distribution 
of Draco with 15 to 20 times higher contrast than 
the previous map of IH95 which has a central den- 
sity contrast of about 4. 

We note that a few small areas (diameter < 1') 
in the central part of Draco show a pronounced 
lack of stars because the source detection and pho- 
tometry has been hampered by the presence of a 
very bright foreground star. These regions have 
either been omitted from the analysis (e.g. in the 
fitting procedures) or else recovered by interpola- 
tion of the stellar density in the surroundings. 

3.1. Surface density maps 

The spatial distribution of the stars in SI and 
S2 is shown in Fig. 4. Individual stellar positions 
are plotted as dots. Overlayed are contour lines 
of constant surface density. The surface density 
has been derived through counts on 3' by 3' grid 
and subsequent weighted averaging within a ra- 
dius of two grid steps in the outer part and one 
grid step in the inner part of the field. The thin 
lines show contours at the level of 1 a above the 
mean background density. Here, a is the rms of 
the background density fluctuations, indicating 
the lowest level at which Draco stars could be- 
come recognizable. This variance corresponds to 
surface densities of Draco stars of 0.035 arcmin"^ 
for SI, and of 0.022 arcmin-^ for S2, i.e. 3 x lO'^ 



of the central surface density of Draco. Any signif- 
icant detection of the Draco population requires 
at least 2 cr above background, marked by the out- 
ermost thick contours in Fig. 4. Further contours 
are drawn at 5 cr and higher in order to reveal the 
shape of the galaxy at different radii. 

Both samples show that the overall distribution 
of Draco stars is approximately ellipsoidal. With 
respect to the innc;r part the outer contours at 
< 5(7 appear to be slightly deformed or shifted to 
the south-west (see §4). Also, the contour lines for 
S2 are in general less smooth than those for SI. 
Nevertheless, there are no clear deviations from 
an elliptic shape down to the 2 a contour. This 
contour reaches out to 30' angular distance from 
the center of Draco. 

Strong departures from an elliptic shape are 

visible only at the 1 a level. Here, the map for the 
sample S2 shows more extended structures than 
the map for SI. There appear to be weak exten- 
sions in four directions out to radii of about 45'. 
The ones pointing to the north-east and south- 
west are filamentary while those in the south-east 
and the north-west are broader and may be de- 
scribed as a bar-like. Beyond 45' distance from the 
center of Draco the field is sprinkled with numer- 
ous patches of 1 cr overdensities. About one third 
of them are visible also on the 2cr level. These 
patches arc distributed rather uniformly, and do 
not form large connected structures. Altogether 
they cover an area of no more than 30% of the to- 
tal field. This suggests that these patches simply 
reflect random fluctuations of the residual field 
star population. A closer look at the statistics of 
the outer fields supports this view. In Fig. 5 we 
show histograms and cumulative distributions of 
the local S2 field star densities measured in non- 
overlapping grid cells of 6' x 6' and 9' x 9' beyond 
45' distance from the center of Draco. The open 
circles show Poisson distributions for the observed 
mean values. The measured fluctuations are in 
close agreement and indeed statistically consis- 
tent with the corresponding Poisson models. A 
Kohnogorov-Smirnow test shows that the opposite 
hypothesis, namely that the empirical distribution 
is different from the Poisson model, has in both 
cases a probability of not more than 2%. Thus 
we conclude that the number of grid cells with 
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surface densities of e.g. 1 cr or 2 cr above the mean 
does not significantly exceed the number expected 
from Poisson noise of a random distribution, (cf. 
Fig. 5b). Hence there is no evidence, and indeed 
httle room, for overdensities that can be ascribed 
to Draco. Another way to show that the density 
enhancements around Draco in tlic maps of Fig. 4 
are random fluctuations is by comparing the ob- 
servations with a random number generated star 
field of equal mean density. Fig. 6 shows such a 
simulated field corresponding to sample S2. The 
surface density has been derived using the same 
pixel grid and smoothing and the density contours 
are drawn on the same levels as for S2 in Fig. 4b. 
The la and 2a contours in the simulated map 
show structures that arc very similar in size and 
number to those in the maps drawn from the ob- 
servational data. 

Density fluctuations of the field star population 
must of course also occur in the region where the 
Draco galaxy is seen. The superposition of these 
fluctuations with the low density of Draco stars in 
the outermost parts of the galaxy is likely to cre- 
ate distortions in the 1 a contour around Draco. 
This is also nicely revealed by the simulated ex- 
ample of Fig. 6. Since the four weak extensions 
which arc visible in Fig. 4b and which have been 
described above are not substantially larger than 
the sizes of the 1 a patches in the outer field it is 
very likely that these extensions are also due to 
fluctuations in the density of residual field stars. 
At least, these structures are no significant indi- 
cation of the presence of Draco stars. 

In summary, we find that down to the level of 
2 a above the background density of our maps the 
Draco dSph has a fairly regular shape, and that 
there are no detectable traces of tail-like exten- 
sions down to 3 X 10^"^ of the central density. 

3.2. Mean surface densities along major 
and minor axis 

The stringent non-detection of tidal tail fea- 
tures on our new Draco maps contradicts the pre- 
vious conclusions of SKH97, who derived mean 
surface densities in fields of 23' x 23' along the 
principal axes of Draco. For a direct comparison 
with the SKH97 results, we determined mean sur- 
face densities from our data in fields of the same 



size and position along the same axes (major axis 
at position angle 82°). These are plotted in Fig. 7 
as a function of position along the major and mi- 
nor axes with respect to the center of the galaxy. 
As in SKH97, the mean density in each field is nor- 
malized to the mean density in the central field, 
and an estimate of the foreground density is sub- 
tracted. Our estimate of the mean foreground 
density, however, relies not only on one remote 
field as in SKH97, but on all stars with more than 
2° distance from the center of Draco. The fields 
at 20' mean angular distance from the center of 
Draco have scaled densities of about 0.15 with re- 
spect to the mean density of the central field, com- 
pared to 0.12 reported by SKH97. In contrast to 
SKH97, we do not observe significant overdensi- 
ties along the major axis far from the center of 
Draco in either of our two samples. In the more 
distant fields, the mean density agrees with the 
density level of foreground stars within the statis- 
tical errors. There is no general diff'erence between 
the mean densities measured along the major axis 
and those measured along the minor axis beyond 
a distance of 40' from the center. However, there 
are of coiirsc statistical fliictuations in the mean 
densities of the individual fields along both axes. 
Fig. 7 shows that the reference foreground field for 
SKH97 (i.e., at 1° north of Draco) is the one with 
the lowest mean density in our series of fields. We 
would obtain significant overdensities in most of 
the other fields if we chose this particular field as 
a reference. This suggests that the overdensities 
reported by SKH97 reflect to a large extent an un- 
derestimated foreground density. The resolution 
of this apparent discrepancy shows once more the 
power of wide area data. 

4. Size and structure of the Draco dSph 
galaxy 

To quantify the size, shape and orientation of 
the Draco dSph we fit two-dimensional models 
to the observed surface densities. Herein we use 
the empirical King profile (King 1962, hererafter 
K62; Eq. 3a), the theoretical King profile (King 
1966, hereafter K66), and a generalized exponen- 
tial (Sersic 1968; Eq.3b): 

S(r) = El • ( ^ ^ - , ^ \ (3a) 
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S(r)=Soexp(-(r/ro)") (3b) 

In all cases r = \J x? + y^/ (1 — e)^ is the el- 
liptical radius and x the coordinate aligned with 
the major axis. The profile Eq. 3a has a finite 
limiting radius rt which is meant to describe tidal 
truncation. It is often used in simple models of 
globular clusters. The theoretical models of K66 
have projected density profiles that resemble those 
of Eq. 3a, but with a somewhat shallower decline 
towards the cut-off radius. These profiles cannot 
be expressed explicitly and were generated by nu- 
merical integration following the description by 
Binney & Tremaine (1987, p.232). The profile 
Eq. 3b is of infinite extent and has been found 
to fit elliptical galaxies in a wide range of abso- 
lute magnitudes (Young & Currie 1994, Jerjen & 
Binggeli 1997). The three types of models were fit 
to the surface density distributions of SI and S2 
by weighted least squares. The surface densities 
were sampled on a 3' x 3' grid of non-overlapping 
cells. A few cells which fell on positions that are 
blended by bright stars were excluded from the 
fit. Weights were derived from Poisson noise, with 
the expected number of stars given by the fitted 
model. The weights were iteratively adjusted to 
yield a consistent solution. The best-fit param- 
eters are given in Table 2 and the best fits are 
shown in Figs. 8 and 9. 

The values of xV(^ - /) in Table 2 as well 
as the residuals plotted in Figs. 8b, 8d and the 
profiles in Fig. 9 demonstrate that the three types 
of models all provide good fits. Moreover, the 
model parameters obtained from the two different 
data samples are consistent within their formal er- 
rors. Hence the results do not depend critically on 
the particular choice of the color-magnitude filter. 
The best-fit ellipticity is 0.29 ± 0.03, confirming 
earlier results of IH95 and Hodge (1964). The 
major axis position angle is 88° (±3) (J2000.0), in 
slight disagreement with the earlier determination 
of IH95 (see Table 2). 

The most interesting parameters are those de- 
scribing Draco's size. Here our results differ sub- 
stantially from those given by IH95. The fit of the 

'(iV — /) denotes the degrees of freedom of the fit and is 
K. 4350. 



K62 model yields for both samples a core radius 
that is about 15% smaller and, more importantly, 
a limiting, or tidal, radius that is about 40% larger 
than the previous standard value of IH95 (see Ta- 
ble 2). The latter has important consequences: (1) 
It shows that Draco is more extended and there- 
fore also more massive than previously assumed 
(see §5). (2) It shows that some of the previous 
detections of Draco stars beyond the former stan- 
dard radius of 28f3 are in fact not cxtratidal if 
we define the (actually unknown) tidal radius of 
the galaxy as usual by the limiting radius of the 
best-fit King profile. 

To demonstrate the clear evidence for a larger 
spatial extent of Draco in our data, we show in 
Fig. 9 the observed radial density profiles for SI 
and S2. The densities are tabulated in Table 3. 
The profiles were derived from star counts in el- 
liptical annuli. The logarithm of the mean density 
above background is plotted versus the logarithm 
of the rms of the outer and inner radius of each an- 
nulus. The radii refer to the major axis. The mean 
background density was obtained as one of the pa- 
rameters in the least-squares fitting of the two- 
dimensional model distributions where it is deter- 
mined by the surface density data in the outer 
parts of the field. Error bars at each bin indicate 
the statistical count rate uncertainty (quadrati- 
cally combined uncertainties of S and E^g) and 
show that all bins except the outermost one have 
significant densities of Draco stars. Note that the 
counts in these bins trace the profile of the galaxy 
over three orders of magnitude and thus reflect a 
ten-fold increase in the dynamic range compared 
to IH95. Three relevant model curves are over- 
layed as lines in Figs. 9a and 9b: two King mod- 
els (dashed lines) and one generalized exponen- 
tial model (solid lines). Clearly, the K62 model 
of IH95 (short dashed lines) fails to describe our 
data, as tlic;re are four consecutive radial bins with 
significantly higher surface density. The new best- 
fit K62 models (long dashed lines) represent the 
observed densities (within their statistical uncer- 
tainties) down to the outermost bin. There are no 
hints for the existence of a break in the slope of the 
radial profile. The models with the K66 surface 
density profile fit the data with the same overall 
quality but provide a better approximation to the 
outermost significant bin (34' < r < 40') since the 
decline of the density profile is less steep. Extra- 
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polation with the best-fit K66 profile thus yields 
a limiting radius of about 50' instead of about 40' 
from K62. The best-fit generalized exponential 
model is characterized by a radial scale length of 
ro = 7f3 (SI) or 7.'6 (S2) and an exponent n = 1.2. 
This profile fits the outermost significant bin in the 
same way as the best-fit K66 models. For sample 
S2 the fit of the generalized exponential even at- 
tains a slightly better overall x^- This reveals that 
there is no compelling evidence for a real cut-off 
in the radial profile of the Draco dSph within the 
current observational limits. In other words, de- 
spite the improved contrast of our density maps 
a distinct edge of the galaxy is possibly still un- 
detected. Our result of rt — 40' from the best-fit 
K62 model may thus underestimate Draco's real 
radius but it puts at least a firm lower limit on 
any possibly existing true density cut-off. 

The model fits provide a useful reference to 
identify peculiarities in the shape of the observed 
distribution of stars. In Figs. 8a and 8c we com- 
pare the contours of the observed density distri- 
bution with the elliptic contours of the best-fit 
exponential model. Apart from the filamentary 
and barlike extensions on the 1 <t level (see § 3) 
one recognizes a less extended box-like deforma- 
tion of the galaxy with respect to the model in 
the south-western quadrant. This feature shows 
up in both samples and is present also in the 2 a 
and 3 a contours. The residual maps (Figs. 8b and 
8d) trace these deformations to two patches of en- 
hanced surface density which are located about 
25' south and southwest of the center of Draco. 
Their 2(t significance might indicate real irregular- 
ities in the stellar distribution of Draco. In order 
to quantify the amount of asymmetry we counted 
the number of stars in two opposite sectors sepa- 
rated by a straight line through the center of the 
galaxy (bisector method). From the (background- 
subtracted) numbers ni, n2 in the two sectors one 
can define an asymmetry parameter A as: 



ni + n2 

For A << 1, the uncertainty of ^ is = 
l/V'^i + "■2- In the inner part of Draco (r < 15') 
both SI and S2 yield no significant asymmetry 
{A < aA, CTA 2.5%) for all position angles of 



the bisector. For 15' < r < 40' we obtain maxi- 
mum values of A = (6.1 ±3.4)% for sample SI and 
A = (9.0 ± 4.6)% for sample 82, both with bisec- 
tor position angles of about 150°. The outer part 
of Draco thus seems to have an asymmetric stel- 
lar distribution at the 10% level {2a significance). 
Note that any tidal deformation is expected to be 
symmetrical with respect to the center. 

5. Mass, luminosity and global mass-to- 
light ratio 

5.1. Estimating the total mass 

The new values for the spatial parameters Tc, rt 
of Draco lead to a revised estimate of the total 
mass M of the Draco dSph. We estimate M in 
the usual way by means of the mass of an equiv- 
alent single-component King sphere. Following 
K66, this mass is given by: 



M=^ni,{f3aiosf (5) 

Herein n denotes the normalized mass (see 
Eq. 38 of K66) which is only a function of the con- 
centration of the system, rc the linear core radius, 
aios the observed line-of-sight velocity dispersion, 
and (3 a correction factor to link the observed dis- 
persion to King's parameter ao. We recall that 
this way of estimating the mass of a stellar system 
relies on three basic assumptions: (1) dynamical 
equilibrium, (2) proportionality of mass and light 
distribution, and (3) the adequacy of a spherically 
symmetric and isotropic model in general and a 
King model in particular. 

Our values for rc and from the best- 
fit K66 profile imply a concentration of c = 
log(rt/rc) = 0.77 ± 0.03. The normalized mass 
thus is = 6.8 ± 0.5 (as opposed to ^ = 2.6 
used in III95). The core radius of the dynamical 
model is appropriately chosen as the geometric 
mean rc = rc,obsV^ — e since this preserves the 
projected area of the system when replacing the 
observed ellipsoidal system by the spherical model. 
To determine the relative change in the mass M 
that is due to changes in rc and it suffices to 
evaluate the product rc ■ n, using rc in angular 
measure. We thus find that M increases by a 
factor of 2.5 when using the revised values of rc 
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Table 2 

Parameter \alues of best-fit ellipsoidal models of Draco 



Sample 


Model 




Sc 


e 


PA 




n 




n 


X'KN - /) 


O -L 




zuu. uoo 


iccyo Q1 7 


u.zo 


oy 




oy. u 






1 niQ 






±0.006 


±0.003 


±0.01 


±3 


±0.2 


±0.8 








SI 


King66 


260? 055 


+57?917 


0.28 


89° 


8.'3 


49.'6 






1.018 














±0.2 


±1.3 








SI 


Sersic 


260? 055 


+57?916 


0.28 


89° 






7f3 


1.2 


1.017 






±0.006 


±0.003 


±0.01 


±3 






±0.1 


±0.1 




82 


King62 


260? 055 


+57? 915 


0.30 


88° 


7.'7 


40.' 1 






0.995 






±0.008 


±0.003 


±0.02 


±3 


±0.2 


±0.9 








S2 


King66 


260? 055 


+57?915 


0.30 


88° 


8.'7 


49.'4 






0.994 














±0.3 


±1.4 








S2 


Sersic 


260? 055 


+57? 915 


0.30 


88° 






7.'6 


1.2 


0.988 






±0.008 


±0.003 


±0.02 


±3 






±0.1 


±0.1 




IH95 


King62 






0.29 


82° 


9f0 


28^3 
















±0.01 


±1 


±0.7 


±2.4 









Note that the meaning of Tc for K62 and K66 is different and depends on the concentration of the model. 



Table 3 

Radial surface density profiles 









SI 








S2 




^out 




E 


a(E) 




E 


a(S) 




arcmin 




arcmin 




arcmin 


0.0 


1.0 


0.71 


10.250 


2.137 


7. 


.181 


1.795 


1.0 


2.0 


1.58 


8.467 


1.121 


5. 


,087 


0.872 


2.0 


3.0 


2.55 


7.219 


0.802 


4. 


,937 


0.666 


3.0 


4.0 


3.54 


7.448 


0.689 


4. 


,937 


0.563 


4.0 


5.0 


4.53 


7.130 


0.594 


4. 


,388 


0.468 


5.0 


6.0 


5.52 


4.740 


0.438 


3. 


,142 


0.358 


6.0 


7.0 


6.52 


5.142 


0.420 


3. 


,452 


0.345 


7.0 


8.0 


7.52 


3.654 


0.329 


2, 


,723 


0.285 


8.0 


10.0 


9.06 


3.181 


0.198 


2, 


,132 


0.163 


10.0 


12.0 


11.05 


2.380 


0.155 


1, 


,642 


0.129 


12.0 


14.0 


13.04 


1.363 


0.108 


0, 


,941 


0.090 


14.0 


16.0 


15.03 


1.099 


0.090 


0, 


,778 


0.076 


16.0 


18.0 


17.03 


0.629 


0.064 


0, 


,442 


0.054 


18.0 


22.0 


20.10 


0.485 


0.037 


0, 


,278 


0.028 


22.0 


28.0 


25.18 


0.300 


0.021 


0, 


180 


0.016 


28.0 


34.0 


31.14 


0.193 


0.015 


0, 


110 


0.012 


34.0 


40.0 


37.12 


0.165 


0.013 


0. 


,088 


0.009 


40.0 


60.0 


50.99 


0.150 


0.006 


0. 


,078 


0.004 



Notes: r,„ = y0.5(rf„ + r2„,) 



The "background" surface densities are: 
Ebg = 0.1511 ± 0.0020 arcmin-2 (SI) 
Ei,g = 0.0760 ± 0.0014 arcmin-^ (S2) 
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and rt instead of the former ones. If the rest of 
Eq. 5 is to be kept unchanged, former estimates 
of M could be updated by just multiplying with 
this factor. However, we consider it worthwhile 
to assess also the values of the other parameters 
involved in Eq. 5. 

The physical value of Vc requires an estimate 
of Draco's heliocentric distance d. Adopting the 
apparent brightness of the red horizontal branch 
measured by Stetson (1979) {Vrhb = 20.07 ±0.03, 
E{B -V) = 0.03), the empirical HB luminosity 
relation of Carney, Storm & Jones (1992) {My = 
0.15 [Fe/H] + 1.01) and a mean metal abundance 
of [Fe/H] = —2.0 (see the compilation by Lehnert 
et al. 1992) one obtains (m - M)o = 19.26 ± 0.20 
and hence d = 71 ± 7 kpc.^^ This yields a lin- 
ear core radius of Vc = 148 ± 15 pc (geometric 
mean). Prom spectroscopic measurements of 91 
Draco members, the most extensive data set of 
radial velocities in Draco that is currently avail- 
able, Armandroff et al. (1995) derived line-of-sight 
velocity dispersions of 10.7 ± 0.9 km s~^ (complete 
sample) and (8.5 ± 0.7) kms~^(by omission of 
three extreme velocities). Olszewski et al. (1996) 
showed that these results are robust against con- 
tamination by orbital motion of binaries. Most 
of the stars of this radial velocity sample are lo- 
cated within an ellipse with semi-major axis of 
1.5rc. By Monte Carlo simulation of a K66 dy- 
namical model with c = 0.77 we find that a cor- 
rection factor P ~ 1.35 applies if the line-of-sight 
velocity dispersion is measured with stars that 
cover the range r < 1.5 fg. Inserting all these 
values into Eq. 5 we obtain estimates of the total 
mass of the Draco dSph of 3.5 ± 0.7 x lO'^M© and 
2.2 ± 0.5 X 1O''M0 depending on which of the two 
velocity dispersions are used. 

As expected, our result of 3.5 x lO^'M© ex- 
ceeds some of the previous estimates derived on 
the basis of the parameters of IH95, but the in- 
crease is less than by the factor of 2.5 mentioned 
above. Hargreaves et al. (1996), for example, ob- 
tained M = 2.6 X lO^M© with the same method 
of King model mass estimation, using rc and rt 

'The uncertainty of 0.2 mag in the distance modulus ac- 
counts for discrepancies between horizontal branch lumi- 
nosity relations derived by different authors and methods. 



from IH95 and a velocity dispersion very near 
to the 10.7 kms~^of Armandroff et al., but with 
a different choice for /3. Mateo (1998) and van 
den Bergh (2000) give the total mass of Draco 
as 2.2 X 10'' M0 referring to similar values for the 
spatial and kinematic parameters, but without 
presenting their calculations in detail. The mass 
estimate of 4.4 x 1O''M0 given in IH95 is substan- 
tially higher than ours and those by other authors. 
This seems confusing but can be traced to their 
adopted dispersion of 13.2 kms~^, which has been 
invalidated by the measurements of Armandroff et 
al. and Hargreaves et al. Our revised estimates of 
Draco's mass benefit from an improved knowl- 
edge of the spatial as well as of the kinematic 
parameters and hence are more precise. Never- 
theless, their validity remains limited to the point 
that the underlying basic assumptions (see second 
paragraph of § 5.1) must hold, at least approxi- 
mately. 

The assumption that "mass follows light" might 
be wrong in the sense that the dark matter com- 
ponent which is likely to exist and to dominate 
the matter content of Draco (see §5.2 and §6.2) 
could have a less concentrated distribution than 
the luminous matter, as is often found to be the 
case in stellar systems that contain dark matter. 
Thus it is interesting to see how the estimate of 
the total mass changes under such circumstances. 
In the most extreme case the dark matter compo- 
nent would have a constant density pdm (out to 
radius rt) and the stars would move as test par- 
ticles in the harmonic potential of this homoge- 
neous dark matter distribution. Here it is con- 
venient to approximate the space density of the 
stars in the core region of Draco by a Gaussian ra- 
dial density profile, i.e., p = po exp(— (r/ro)^/2), 
with scale length ro = rc/j and 7 ~ 1.42 for 
c 0.77. From the Jeans equation (see Binney 
& Tremaine 1987, p. 204) it follows that a system 
with Gaussian density profile in a harmonic po- 
tential has a velocity distribution with constant 
dispersion a = \J AirGpoM /'^ (see alsO Pryor & 
Kormendy 1990). Setting a = aios, the mass M of 
the dark component out to radius rt and the ratio 
between this M and the King mass M of Eq.5 are 

^=i4^L (6) 
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M^4n fjV (nV 1 

M 9 \0J \rcj M ^ ^ 

Putting the appropriate parameter values into 
Eq. 7 one finds M/M ~ 45. This example demon- 
strates that the mass estimate derived by Eq. 5, as- 
suming that all matter in Draco follows the distri- 
bution of the stars, should be regarded as a lower 
limit because a more shallow distribution of the 
dark matter component would result in a consid- 
erably higher total mass. 

5.2. Luminosity and mass-to-light ratio 

We determined the total luminosity of the 
Draco dSph by (1) summing up the i band fluxes 
of all possible Draco members in the magnitude 
interval 16.0 < i* < 21.5, (2) subtracting the ex- 
pected contribution of flux from foreground stars, 
and (3) adding an estimate of the missing flux from 
Draco stars at i* > 21.5. The flux integral was 
calculated in the region of the color-magnitude 
plane that is outlined by the solid line in Fig. 10, 
and for major axis radii varying from 10' to 50'. 
The foreground contribution was estimated by an 
analogous integration in a ring around Draco ex- 
tending from 1° to 2°. The foreground-corrected 
luminosity function of Draco is shown in Fig. 10b. 
After foreground subtraction the integrated mag- 
nitude for Draco converges to rrii = 11.09 ± 0.03 
at a radius of r = 20'. At r > 20' the integrated 
magnitude deviates by not more than ±0.03 mag 
from this value. The uncertainty of mj was esti- 
mated by splitting the fields into four quadrants 
and calculating the flux integral for each quadrant 
separately. The resulting rms was 15 %, implying 
an uncertainty of 0.08 mag for the determination 
of rui from the entire field. As a further check we 
mention that an erroneous inclusion or omission of 
up to three stars at the bright end of Draco's giant 
branch (i* = 16) would change rrii by not more 
than 0.04 mag, which is only half of the quoted 
error bar. 

The missing flux from faint stars beyond our 
magnitude cut-off was estimated by analysing 
deep HST observations of a small field near the 
center of Draco (Grillmair et al. 1998). We re- 
trieved these observations from the HST archive 



and reduced them with the software package HST- 
phot (Dolphin 2000). The resulting photometry 
covers the magnitude range 18.5 < I < 26.0 and 
has an overlap of 2.8 mag with the SDSS data in 
the range 18.5 < / < 21.3 or 18.9 < i* < 21.7. 
The luminosity function derived from the HST 
data is shown by the dashed line in Fig. 10b. We 
integrated the / band counts from HST between 
7 = 21.3 and I = 26.0 and extrapolated the inte- 
gral to / = 27.75 assuming a constant slope of 
the luminosity function. We then scaled the re- 
sult with the ratio between the number of Draco 
stars in the (global) SDSS sample and the num- 
ber in the (local) HST sample, both counted in 
the overlap interval. This yields an estimate for 
the missing flux in the SDSS data of 0.59 (±0.08) 
mag. Correcting for the faint stars in Draco thus 
leads to rrii = 10.49 ± 0.10 for the total mag- 
nitude. By subtraction of the distance modulus 
((m — M)o = 19.26, see § 5.1) nii translates into 
a total absolute magnitude of Mj = —8.77 ± 0.20. 
The uncertainty of the absolute magnitude is dom- 
inated by the uncertainty of the distance modulus 
of 0.2 mag. 

For conversion to solar units we used basic 
photometric data for the Sun given in Landolt- 
Bornstein (1981) and an empirical transformation 
relation between standard photometry and SDSS 
photometry (Grebel et al. 2001): 

i* = 0.33(±0.02) + Ic + 0.19(±0.02) {V - Ic) (8) 

This relation was derived by comparing SDSS 
data with Vic photometry in the photometric 
standard fields of Stetson (2000). For the Sun we 
have Mv,Q = 4.87 and a color index iV — /)© = 
0.81 in the Johnson system, which, according to 
the transformation formulae of Fernie (1983), cor- 
responds to iy — Ic)q = 0.64 in the Johnson- 
Cousins system. Putting this into Eq. 8 one ob- 
tains Mi*^Q = 4.68. Our result for the luminosity 
of Draco thus reads {L/LQ)i = (2.4 ± 0.5) x 10^. 
This is in good agreement with earlier determi- 
nations of (L/Lq) made in a different wavelength 
range, e.g. by IH95 who give the luminosity of 
Draco as {L/Lq)v = (1.8±0.8) x 10^ We empha- 

Deep HST observations of globular clusters by Piotto et al. 
(1997) show that the maximum of the luminosity function 
normally lies at about Mj = 8.5, i.e. / = 27.75 for Draco. 
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size that the determination of the hiniinosity does 
not depend on critical assumptions and models 
and is in this sense more general than the deter- 
mination of the total mass. 

Combining the luminosity and the mass esti- 
mates of § 5.1 we find total mass-to- light ratios 
M/Li of either 146 ± 42 or 92 ± 28 in solar units. 
For comparison, the values of M/ Ly given by Ar- 
mandroff et al. (1995), were 90 and 57 in solar 
units (with respect to the V band). As for the 
values for the total mass, the mass-to-light ratios 
are valid only on the proviso that the assumptions 
mentioned in § 5.1 hold. Note however, that the 
true mass of the system would need to be smaller 
by a factor of 30 or more in order to achieve a 
normal mass-to-light ratio of 1 - 3 in solar units. 
On the other hand it has been shown that a dis- 
tribution of dark matter that is less concentrated 
than the distribution of the stars would result in 
an increase of the total mass by up to a factor 
of 45 and thus lead to a further increase of the 
mass-to-light ratio. 

An additional argument for a high mass-to- light 

ratio comes from the size limit that is imposed by 
the galactic tidal forces. This argument is inde- 
pendent from the above direct estimate of Draco's 
total mass since it relies only on the observed max- 
imum spatial extent and the luminosity. For a 
satellite of mass M orbiting in a spherically sym- 
metric potential at (constant) galactocentric dis- 
tance R there exists a well-defined theoretical size 
limit that is given by the Lagrange radius char- 
acterizing the force balance between internal grav- 
ity, differential centrifugal forces and the galactic 
tidal field. Taking a logarithmic potential with cir- 
cular velocity Vc = 220 kms^^as a standard model 
for the gravitational potential in the Galactic halo, 
the Lagrange radius is given by the relation: 



Fig. 11 shows rx, as a function of galactocen- 
tric distance R for different satellite masses M 
(dashed lines) and different mass-to-light ratios 
(solid lines), adopting the observed luminosity 
{L/LQ)i = 2.4 X 10^ of Draco. The fact that 
no signs of tidal perturbation have been detected 



suggests that at least the currently observed part 
of Draco's stellar distribution does not reach be- 
yond the Lagrange radius. Therefore we have the 
condition rt < r^. The measured limiting radius 
Tt is drawn as the hatched bar in Fig. 11. The 
galactocentric radial velocity of Draco which 
is near —100 kms^-'^, indicates that Draco is not 
on a circular orbit. Nevertheless, one can use tl 
for obtaining an approximate size limit. The ap- 
propriate distance R at which to take is in this 
case however not evident; it might be the pericen- 
tric distance Rp rather than the current distance 
i?o • Since Rp cannot be estimated without reliable 
information on the tangential velocity, one has to 
consider a wide range of distances R< Rq. Fig. 11 
reveals that for any such choice of i? a normal 
mass-to-light ratio of M/L < 3 yields Lagrange 
radii that are substantially smaller than the ob- 
served limiting radius rt for which we here take the 
conservative value of 40' or 820 pc. Hence a nor- 
mal value of M/L is ruled out within the context 
of our assumptions. Possible values of M/L in the 
sense that rt < tl holds for distances well below 
Ro lie near or above M/L = 30. With M/L « 100, 
i.e. a satellite mass of about 3 x 10^ Mq as derived 
by direct mass estimation with Eq. 5, the observed 
size of Draco remains below the theoretical limit 
of down to galactocentric distances of 20 to 
25 kpc. This means that a mass-to-light ratio of 
about 100 or higher can prevent tidal stripping of 
stars from the outermost part of Draco even in rel- 
atively close approaches to the Milky Way. It must 
however be noted that this is a stationary picture. 
In reality a passage of the Galactic center on an 
elongated orbit could have the additional effect of 
heating the stellar population of the satellite due 
to the non-stationarity of the tidal forces. This 
kind of influence on the internal kinematics of the 
satellite system is not accounted for by the simple 
Lagrange radius criterion. 

6. Discussion and conclusions 

6.1. The missing extratidal component 

From the results of previous studies, in partic- 
ular those of IH95 and SKH97, one might have 
expected that a high contrast map of the Draco 



Note that the line of sight from the Sun to Draco is close 
to Draco's galactocentric radius vector 
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galaxy down to magnitude 22 would show clear 
signs of tidal features. Our analysis of the SDSS 
data shows that this is not the case. Structures 
which could be interpreted as indicating tail-like 
extensions appear only on the 1 a level of our most 
sensitive map, i.e., at 2.8 x 10~^ of the galaxy's 
central density. Since these structures go into four 
different directions their interpretation in terms of 
tidal tails is problematic. The more bar-like struc- 
tures seen towards the south-east and the north- 
west would appear as the preferred candidates for 
such tails. However, by statistical analysis of the 
density variations in the surrounding field we find 
that the observed amount of structure at 1 a above 
background density is easily explicable by normal 
Poisson noise. Therefore, none of the extensions in 
the vicinity of the main body of Draco nor any of 
the overdense patches seen at larger angular dis- 
tance are likely to be an external part of its stellar 
population. 

On the other hand, the radial profile of the 
surface density shows that the ordinary spatial 
extent of the Draco galaxy is substantially larger 
than was revealed by previous measurements with 
lower contrast. Going down to about 10~^ of the 
central density, the observed radial density pro- 
files directly reveal that Draco extends to about 
40'. The same is found by fitting ellipsoidal mod- 
els to the observed surface density distribution. A 
K62 profile yields a best-fit rj of 39.'6 (SI) and 
40.'1 (S2), which is 1.4 times the previous deter- 
mination (all radii referring to the major axis). 
Using the profile of the K66 theoretical models 
one obtains an even larger estimate of Vt = 49 f5. 
The observations do however not provide particu- 
lar support for a truncated radial profile. We find 
that a model with generalized exponential profile 
and exponent n ~ 1.2 fits the surface density dis- 
tribution equally well or even slightly better. Thus 
it is possible that the density profile of the galaxy 
continues further outwards and that its true limit 
lies at densities that are below the current obser- 
vational noise level. 

In the light of these results, Piatek et al.'s re- 
cent evidence for Draco stars immediately beyond 
the border of the standard ellipse of semi-major 
axis 28;3 (fields El, Wl, SEl of Piatek et al. 2001) 
has to be understood as a detection of the ordi- 
nary Draco population and is not an indication 



of an extratidal component. Piatek et al. claim 
that there is also weak evidence for Draco stars in 
their fields E2 and W2, i.e. at angular distances of 
about 1.3° from the center. However they say that 
one may disagree with this view because their data 
are such that they can neither prove nor disprove 
the presence of Draco stars at these larger radii. 
Our results indeed suggest that there is no Draco 
component at large radii down to surface densities 
of lO^"* of the central density. In the elliptical an- 
nulus with semi-major axis radii between 40' and 
60' we find mean surface densities above back- 
ground that are < 6 x 10"'' (SI) and < 1 x 10"^ 
(S2) of the central density and compatible with 
zero within the statistical error. IH95 however 
reported densities above background of the order 
of 10~^ times the central density in the distance 
range from 40' to 70' (cf. panel 2 of their Fig. 2). 
This result has often been regarded as a major 
piece of evidence for the existence of an extratidal 
Draco component. With a limiting magnitude of 
i? ~ 21 the IH95 study of Draco does not com- 
prise fainter stars than our SDSS samples. On 
the other hand, above the limit of i* = 21.7 our 
sample SI is unlikely to miss a substantial frac- 
tion of the Draco population since the filter mask 
is sufficiently wide (see Fig. 3c) to include all rel- 
evant parts of the color-magnitude diagram. Yet, 
the color-magnitude filter provides us with the ad- 
vantage of a much higher contrast. Our counts of 
specifically "Draco-like" stars thus prove that sur- 
face densities as high as those reported by IH95, 
which are ten times higher than in our samples, 
cannot be due to Draco members. An extratidal 
Draco component with a surface density as given 
by IH95 is clearly ruled out. Instead, we attribute 
the counts of IH95 beyond 40' to residual fore- 
ground field stars. 

The results of the study by SKH97, which 
also favored an extratidal Draco component, have 
turned out to be at least doubtful. A detailed 
comparison, for which mean densities were deter- 
mined in the same fields as observed by SKH97, 
revealed that their detection of overdensities may 
be due to an accidental underestimation of the 
density of foreground stars and thus spurious. We 
find that the field that was used to set the fore- 
ground reference in SKH97 has an extraordinarily 
low surface density of "Draco-like" stars which 
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is not representative of other and more distant 
fields. With respect to the mean density of stars 
at more than 2° distance from Draco we do not 
find significant overdensities in the SKH97 fields. 

In total, there remains no convincing evidence 
for an extratidal Draco component. This however 
does not necessarily mean that such a component 
is completely absent. For the magnitude range 
down to i* = 21.7, i.e. the giant branch, horizon- 
tal branch and subgiant branch, our study yields 
upper density limits of 3 x 10~^ times the central 
density (from the two-dimensional distributions) 
and of 10~^ times the central density (from the 
radial profiles). The latter limit corresponds to a 
surface brightness of about 33.5 magarcsec^^. In 
a recent numerical study of the dynamical evolu- 
tion of dwarf galaxies in a Milky Way halo poten- 
tial Mayer et al. (2001) found tidal streams with a 
maximum surface brightness of 30 magarcsec"^. 
An extratidal Draco component on this level of 
surface brightness would have shown up in our 
data and can therefore be ruled out. Extratidal 
extensions with surface brightness below the above 
given limits may however exist. N-body simula- 
tions of the tidal disruption of dwarf galaxies by 
Oh, Lin & Aarseth (1995) have indeed produced 
tidal tails with surface densities between 10~^ and 
10~^ of the central density of the dwarf. Thus it 
may be that despite the significant improvement 
in surface density contrast our maps of Draco 
are not yet sensitive enough to reach the level 
at which tidally induced extensions turn up. At 
magnitudes beyond our cut-off, nothing is known 
about signs of tidal perturbation because no ob- 
servations are available so far except for the very 
center of Draco. Since one docs not expect mass 
segregation to be present in such a low-density 
system like Draco tidal effects should in this case 
not depend on the brightness of the tracer popu- 
lation. Therefore the above given limit in relative 
surface density should hold also for fainter stars. 
However, fainter stars, i.e. in the region of the 
main-sequence turn-off and below, have the ad- 
vantage of being much more numerous than those 
on the giant and subgiant branch (sec the lumi- 
nosity function Fig. 10b). On the other hand the 
number of field stars with colors in the range of 
the Draco turn-off will not increase accordingly, 
because the contaminating stars in this part of 



the CMD will be halo turn-off and main sequence 
stars at about the same distance as Draco. Thus, 
from the point of number statistics and contrast, 
Draco stars at and beyond i* = 23 mag appear to 
be much better suited for an eventual detection 
of an extratidal component. The goal of future 
wide-field studies of Draco must therefore be to 
push the magnitude limit down to below the main 
sequence turn-off. 

Since there is no evidence for an extratidal com- 
ponent around Draco (down to the above limit 
of surface density) the question arises why other 
Milky Way dSph's such as Sagittarius, Ursa Ma- 
jor and Carina do show tidal effects (see §1) while 
Draco does not. The clearest evidence of tidal 
perturbation undoubtedly exists for the Sagittar- 
ius dwarf. The case of Sagittarius is however not 
comparable to any of the other known Milky Way 
companions since it is much closer to the Galactic 
center {R ~ 19 kpc for Sagittarius, i? > 66 kpc 
for other M\\ clSph's). Thus for Sagittarius one 
can directly infer from its current position that 
it must be subject to strong tidal forces whereas 
this is not certain for other more distant dwarf 
spheroidals as long as their orbits are unknown. 
In other words, it is by no means clear whether 
Draco or any of the other dSph's has penetrated 
as deeply into the Galactic potential during its 
history as we presently observe it for the Sagit- 
tarius. The Ursa Major dSph is comparable to 
Draco in that its present galactocentric distance 
{R ~ 66 kpc, Mateo 1998) is similar to that of 
Draco, that it has a very low luminosity, and that 
a very high mass-to-light ratio [M/L ~ 80, Mateo 
1998) is obtained via conventional mass estima- 
tion based on the central velocity dispersion. On 
the other hand, the morphology of the Ursa Major 
dSph is quite different from that of Draco because 
it is much more flattened (e = 0.56, IH95) and 
reveals a clumpy internal structure (1H95; Kleyna 
et al. 1998). These peculiarities suggest that the 
Ursa Major dSph has undergone much stronger 
tidal interactions with the Galaxy than Draco. It 
is then not so surprising that Ursa Major (in con- 
trast to Draco) has a supposedly extratidal halo 
with a surface density of about 10~^ of the central 
surface density as found by Martinez-Delgado et 
al. (2001). The reason might be that Ursa Major 
is on a galactic orbit that has smaller pericentric 
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distance than the orbit of Draco. In this case it 
would however be unUkely that Draco and Ursa 
Major belong to the same great circle stream as 
proposed by Lyndcn-BcU & Lyndcn-Bell (1995). 
The key to answering these questions lies in precise 
measurements of the tangential velocities (see end 
of § 6.2). The case of the Carina dwarf spheroidal 
is difficult to understand as long as no information 
on its orbit is available. With a galacto centric dis- 
tance of about 100 kpc (Mateo 1998) it is presently 
much more distant than Draco, so strong tidal per- 
turbations seem a priori less likely. On the other 
hand, Majcwski ct al. (2000) have identified an 
extended distribution of Carina-like giants with a 
surface density on the level of about 10~^ times 
the central surface density and a shallow radial 
profile which is thought to be rather typical of an 
extratidal component. The fact that Carina's to- 
tal mass as estimated in the usual way from its 
velocity dispersion and its surface density profile 
(M/Mq ~ 1.3 X 10'^, Mateo 1998) is a factor of 
two to three lower than the likewise estimated to- 
tal mass of Draco is not sufficient to explain the 
apparently very different amounts of tidal mass 
loss from these objects. One may wonder whether 
the remarkably different star formation histories 
of Carina and Draco perhaps indicates fundamen- 
tal differences in the evolution of the two galaxies 
that could have an impact also on the existence of 
an extratidal component. 

6.2. Draco's structure and basic parame- 
ters 

Apart from setting a new and stringent upper 
limit on the surface density of extratidal stars 
around Draco, the SDSS data provide us with im- 
proved knowledge of the structure and the basic 
parameters of Draco's main body. We find that 
its spatial extent and projected shape is essen- 
tially that of an ellipsoid with a core radius of 
about 160 pc (referring to a decline to about 26% 
of the central surface density), an outer radius 
of at least 820 pc (both along major axis), and 
a radius-independent ellipticity of 0.29. There is 
possibly a 10% deviation from this simple sym- 
metric structure in the outer part. The simplicity 
of the spatial structure of the Draco dSph is a 
very important fact, as is the lack of evidence 
for an extratidal component around it, because 
both findings lend support to the fundamental 



assumption that the system is in dynamical equi- 
librium or at least close to it. This is one of the 
key assumptions generally made in order to es- 
timate the dynamical masses of systems like the 
dwarf spheroidals. The shape of Draco's radial 
density profile further suggests that a King model 
should indeed be an acceptable approximation 
of the system. Under these assumptions and in 
combination with kinematic data the spatial char- 
acteristics of the system imply that it has a total 
mass of 2.2 - 3.5x10'^ Mq. If the distribution of 
mass in Draco is less spatially concentrated than 
is indicated by its luminous matter, the total mass 
out to rt could be by up to a factor of 45 higher. 
To check whether the equilibrium hypothesis and 
the dynamical model is really adequate one would 
need detailed information on the velocity profile 
of Draco, which is currently not known. A more 
definitive answer on the question of the total mass 
of Draco must therefore still await future detailed 
kinematic measurements. 

The luminosity of Draco has been determined 
in a straightforward and model-independent way 
by integrating the observed flux on Draco's giant 
branch, horizontal branch, and subgiant branch 
over the complete volume of the galaxy. For the 
residual luminosity from Draco stars beyond the 
magnitude limit of the SDSS sample we also ob- 
tained a direct and almost assumption-free esti- 
mate using supplementary Draco data from HST. 
The resulting total magnitude of rrii* = 10.50 cov- 
ers the complete range from the brightest giants 
at i* = 16.0 down to the presumed maximum of 
the luminosity function at i* = 28.1. Thus there 
remains little room for significant deficits due to 
still unaccounted contributions of light. The total 
i band luminosity is {L/Lq)^ = 2.4 x 10^, clearly 
demonstrating that the Draco dSph is substan- 
tially less luminous than the brightest globular 
clusters in our Galaxy although it is at least a 
factor of 10 more spatially extended (in the sense 
of its linear diameter) than those clusters. 

The low luminosity is in obvious contrast with 
the estimated total mass of Draco. A total mass of 
2.2 - 3.5x1O''M0 yields a global mass-to-light ra- 
tio of 92 - 146 {M/Li)Q with respect to the i band, 
which lies at least 1.5 orders of magnitude above 
the normal values of the mass-to-light in systems 
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with similar content of luminous matter, namely 
globular clusters (1-3 {M/L)q). In the compila- 
tion of Mateo (1998, Table 4) the only other Local 
Group dwarf galaxy with a similarly high Ad / L 
is Ursa Minor, while most other dSphs have esti- 
mated M/L ratios < 40. Taken at face value, the 
above mass-to-light ratio necessarily implies that 
Draco must be dominated by some sort of dark 
matter component. Even if our mass estimate 
were wrong by a factor of 2 - 5, the conclusion 
that a large fraction of the mass in Draco must 
come from a dark component is inescapable. A 
deviation from the assumption that "mass follows 
light", in the sense that the dark matter compo- 
nent in Draco could - as in other galaxies - be 
less concentrated than the visible matter, would 
further increase the mass-to-light ratio. 

Hence there exist essentially two possibilities: 
(a) The conventional mass estimates must deviate 
strongly from the true mass so that the real mass- 
to-light ratio is normal. This could be the case 
if the system is very far from equilibrium or un- 
bound, (b) The Draco dwarf must contain a large 
fraction of its mass as dark matter, whatever the 
exact quantity of this component may be. The 
fact that the spatial structure of the Draco popu- 
lation does not reveal any sign of a perturbation 
of its internal dynamics makes the first possibility 
unlikely and points towards the second one. 

Using a simple logarithmic halo potential we 
showed that constraints from the Galactic tidal 
field strongly support this conclusion. With a nor- 
mal mass-to-light ratio, which means a total mass 
of not more than 7 x 1{]''Mq, the observed size 
of Draco would clearly exceed the Lagrange limit 
for any reasonable value of the galactocentric dis- 
tance. Hence the outer part of Draco should be 
heavily affected by tidal stripping, which is not 
observed. How much higher the mass-to-light ra- 
tio has to be depends on how closely Draco's orbit 
approaches the Milky Way. Taking twice the value 
of the radial velocity as a guess for the unknown 
tangential velocity yields a hypothetical pericen- 
tric distance of 48 kpc. To reach consistency be- 
tween the measured radius rt and the Lagrange 
limit down to this pericentric distance a mass-to- 
light ratio of M/L k. 30 is required. Values near 
M/L =100 would enable Draco to maintain its 
current size even at pericentric distances of 25 kpc. 



Whether Draco is likely to come that close to the 
Galaxy can only be found out with very precise 
measurements of Draco's absolute proper motion. 
The tangential velocity would in this case be about 
120 kms~^and the (galactocentric) proper motion 
about 0.4 mas/a. The future astrometric satellite 
missions SIM and GAIA are expected to measure 
this motion with high accuracy so that firm con- 
straints on Draco's orbit and hence on the theo- 
retical tidal boundary can then be obtained. 

This paper has focussed on analysing the prop- 
erties of the Draco dwarf spheroidal galaxy as a 
whole. Apart from this the SDSS data also en- 
able us to investigate the stellar populations of 
Draco (i.e., red giants, AGB stars, subgiants, red 
and blue horizontal branch stars, and RR Lyrae 
variables) individually and to compare their dis- 
tribution. This study is currently under way and 
will be described in a future paper (Odenkirchen 
et al. 2001, in preparation). 
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(a- 360.0) COS £ [deg] 



Fig. 1. — Spatial distribution of SDSS point sources (dots) and interstellar extinction (contours) in the region 
of the Draco dwarf spheroidal galaxy. The concentration of dots in the center of the field reveals the dwarf 
spheroidal. In order to avoid overcrowding only 1/10 (!) of the total sample of point sources is plotted. The 
contours indicate extinction levels of 0.08, 0.12, 0.16, and 0.20 mag in Ag* in the order of increasing line 
thickness. Extinctions in u, r, i, and z are obtained by scaling the values of Ag with 1.36, 0.73, 0.55, and 
0.39, respectively. The extinction values are based on the reddening maps of Schlegel et al. (1998). 
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Fig. 2. — Color magnitude diagrams of SDSS stars in the core region of the Draco dSph (elhpse with semi- 
major axis of 9' and elhpticity 0.3 centered on {a, 6) = (260.056°, +57.915°)) using different color indices. 
The diagram of i* versus Ci (panel (c)) optimizes the distinction of the different types of Draco stars, i.e., 
giant branch, asymptotic giant branch, subgiant branch, red horizontal branch, RR Lyrae stars and blue 
horizontal branch. In the "orthogonal" color index C2 (panel (d)) Draco stars of all magnitudes scatter 
around zero. For the definition of the indices ci and C2, see Eq.l. 
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c, [mag] c, [mag] c^ [mag] 

Fig. 3. — Contour plots showing the distribution of stars in the eolor-magitude plane of (ci, i*). (a): Density 
distribution /d of Draco stars from counts in an ellipse of semi-major axis 30', ellipticity 0.3, and position 
angle 90°, centered on the Draco dSph (contribution from field stars subtracted), (b): Density distribution 
/f of field stars from counts at > 2° angular distance from the center of Draco, (c): Lines of constant 
number ratio s = fo/ fp- In all plots the contour levels are equidistant and increase with line thickness. 
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(n- 260.0) cos 6 [deg] 

Fig. 4. — Spatial distribution of photometrically selected Draco candidates. Dots represent individual stars, 
lines are contours of equal stellar surface density. Th(2Siin line shows the level of la above mean background 
density, bold lines show levels of 2a, 5a and higher {a = rms background variation), (a): Sample SI, i.e., 
stars lying within the outermost contour line of Fig. 3c. Contour levels are 0.185, 0.22, 0.31, 0.6, and 1.6 
arcmin"^. The mean density of the background is 0.154 arcmin"^. (b): Sample S2, i.e., stars lying within 
the middle contour of Fig. 2c. Contour levels are 0.097, 0.12, 0.18, 0.3, and 0.8 arcmin"^. The mean density 
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Fig. 5. — Statistics of the local surface density of field stars in the sample S2. (a) Normalized histograms 
of density values measured in pixels of 6' x 6' (upper solid curve) and 9' x 9' (lower solid curve) at angular 
distances of more than 45' from the center of Draco. For comparison, Poisson distributions for the observed 
mean values are drawn (dashed lines) . (b) Cumulative distributions corresponding to the histograms shown 
in (a). The fraction of pixels with densities above a given threshold is seen to agree with predictions for 
Poissonian fluctuations. 
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Fig. 6. — Monte-Carlo simulation of a star field with the same basic parameters as sample S2. The back- 
ground has constant mean surface density of 0.076 arcmin"^, the Draco galaxy is modelled as an ellisoid with 
a generalized exponential profile (see §4 and Table 2) . The surface density has been sampled and smoothed 
in exactly the same way as in the case of the observed samples shown in Fig. 4. The contour level are the 
same as those in panel (b) of Fig. 4. 
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Fig. 7. — Relative mean surface density in subfields of 23' x 23' placed along the conventional principal axes 
of the Draco galaxy. The densities are scaled to the mean density above background in the central subfield. 
The mean background density, determined from stars at more than 2° distance from Draco, is subtracted. 
Panel (a) is for sample SI, panel (b) for sample S2. Open circles show the densities for the fields along the 
major axis (position angle 82°), filled circles are for the fields along the minor axis. The angles along these 
axes increase from south to north and from cast to west, their zero point is in the center of Draco. These 
plots are designed such that they allow direct comparison with Fig. 2 of Smith et al. (1997, SKM97). The 
field that was used as the foreground reference in SKM97 (1° north of Draco) is marked by an open square 
(for details see text). 
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Fig. 8. Contour plots of the observed distribution of stars of Draco, the corresponding best-fit elfiptical 
model distribution, and the residuals of this fit. (a) Surface density contours for sample SI: Thin solid lines 
show the observed surface density (the same as in Fig. 4). The three outermost contours are drawn at levels 
of 1 (T, 2(7 and 5cr above background, a being the rms of the background variations. Thick dot-dashed lines 
show the contours of the best-fit exponential model (Sersic profile with n = 1.2) at the same levels, (b) 
Residuals (O — C) of the fit shown in panel (a), rescaled to the level of the mean background counts by 
means of the square root of the local amplitude of the model. Solid lines show positive residuals, dashed 
lines show negative ones. The contour levels correspond to ±la and ±2a. (c): Same as (a), but for sample 
S2. (d): Same as (b), but for S2. For details on the samples and the models, see text. 
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Fig. 9. — Radial profiles of the surface density of stars in Draco. Data points with errorbars show the profile 
of the observed distribution (i.e., star counts in elliptical rings around the center of Draco, mean background 
subtracted). The lines show the profiles of different models, i.e. a K62 model with the parameter values 
given by IH95 (dotted), the new best-fit King models, K62 (short dashed) and K66 (long-dashed), and the 
best-fit Sersic model (solid, see also Fig. 7). Panel (a) is for sample SI, panel (b) for sample S2. 
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Fig. 10. — The luminosity lunction of Draco from SDSS and HST. Panel (a) specifies the region in SDSS 
color-magnitude space used to construct the bright part of the luminosity function. For comparison the dots 
show the stars in the core region of Draco, as in Fig. Ic. Panel (b) shows the luminosity function (solid 
line) that results from counting the SDSS stars in the region outlined in (a) within the tidal radius of Draco. 
The contribution from field stars has been subtracted. Error bars indicate combined ^Jn uncertainties. The 
dashed line in panel (b) shows the luminosity function of fainter Draco stars derived from 7-band HST 
observations (see text). The magnitude bins arc shifted by 0.4 mag to match the mean offset between ?'* 
and /, the counts are rescaled such that the integrated numbers in the overlap interval (18.9 < i* < 21.7) 
coincide. 
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Fig. 11. — Lagrange radius rj, of the Draco dSph in a logarithmic spherical Galactic potential with circular 

velocity 220 kms"-'^, plotted as a function of galactocentric distance R and for different satellite masses. The 
solid lines assume the observed luminosity of L = 2.4 x lO^L© and mass-to-light ratios of 3, 30, and 100. 
The shaded region indicates the tidal radius of Draco as determined by the fit of the King (1962) profile (i.e. 
rt = 40', corresponding to 820 pc). The current galactocentric distance of Draco is i? = 71 kpc. Its pericentric 
distance Rp must clearly be smaller, since Draco has a galactocentric radial velocity of Vr = —100 kms~^. 
Assuming, e.g., Vt = 200 kms~^for the unknown tangential velocity component, the pericenter would be at 
Rp = 48 kpc. 
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